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ABSTRACT 

We present the ultraviolet-optical-infrared spectral energy distribution of 
the low inclination novalike cataclysmic variable V592 Cassiopeiae, including 
new mid-infrared observations from 3.5-24 yum obtained with the Spitzer Space 
Telescope. At wavelengths shortward of 8 /im, the spectral energy distribution of 
V592 Cas is dominated by the steady state accretion disk, but there is flux density 
in excess of the summed stellar components and accretion disk at longer wave- 
lengths. Reproducing the observed spectral energy distribution from ultraviolet 
to mid-infrared wavelengths can be accomplished by including a circumbinary 
disk composed of cool dust, with a maximum inner edge temperature of ~ 500 
K. The total mass of circumbinary dust in V592 Cas (~ 10^^ g) is similar to 
that found from recent studies of infrared excess in magnetic CVs, and is too 
small to have a significant effect on the long-term secular evolution of the cat- 
aclysmic variable. The existence of circumbinary dust in V592 Cas is possibly 
linked to the presence of a wind outflow in this system, which can provide the 
necessary raw materials to replenish the circumbinary disk on relatively short 
timescales, and/or could be a remnant from the common envelope phase early in 
the formation history of the system. 
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Introduction 



V592 Cassiopeiae was discovered by lGreenstein et al.l (llQTOl ). Its blue color and spectral 
characteristics identified it as a disk-accreting novalike cataclysmic variable (CV). Its optical 
spectrum displays strong He II 4686 A and C III/ N III Bowen blend e mission, with weak 
Balmer and He I emission lines ( iDownes et al.lll995l ). iHuber et al.l (119981 ) performed the first 
time-resolved optical spectroscopic study of V592 Cas and used radial velocity measurements 
to determine an init ial estimate of the o rbit al period, Pnrh = 2-47 h r. This was later refined 
to Porb = 2.76 hr bv lXavlor et~aD Jl998h and Iwitherick et al.l J2003h . Although photometric 



obser vations showed that it does not eclipse (orbital inclination of i = 28°; 



1998 ). V592 Cas does show both positive and negative permanent superhumps (ITaylor et al. 



Huber et al. 



19981 ). Optical and ultraviolet (UV) spectroscopic observations demonstrate d the presence of 



a wind, which maiiifests primarily as a P Cy g ni profile in the emis sion lines (IWitherick et al. 



20031 : IPrinia et al.l l2004l : iKafka et al.ll2008f ). iKafka et al.l (120081 ) showed that the wind is 



variable and reaches velocities up to f„ind ~ 5000 km s ^. 

The secular evolution of interacting binary stars is driven by the transfer of matter 
between the stellar components and out of the system, and the corresponding redistribution 
and loss of angular momentum. However, the observed orbital period distribution of the 
CV population (which provides a snapshot of the evolutionary process) shows several dis- 
crepancies compared to the standard two-mechanism (gravitational radiation and magnetic 
braking from the secondary star) model for angular momentum loss. For example, pop- 
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shorter (~ 65-70 min; e.g., 
the observed minimum of ~ 
two mechanism model predicts a tight relationship between orbital period and mass trans- 
fer rate, whereas the CV population shows a s pread of up to several orders of m agnitude 
in mass transfer rate at a given orbital period (jPattersoru 1 19841 : iKolb et al.ll200ll ). An ad- 
ditional angular momentum loss mechanism, which would cause CVs to reach a minimum 
orbital separation (and, hence, period) sooner during their secular evolution and contribute 
to the spread in mass transfer rates, could account for this. 

In this paper, we use new infrared (IR) observations of V592 Cas obtained from the 
Spitzer Space Telescope to explore a mechanism that has been proposed as an additional 
source of angular momentum loss in CVs: gravitational torques on the inner binary from 



3 



a circumbinary disk , likely composed of cool dust (ISpruit fc Taaml 1200 ll : iTaam et al.l 12003 



Willems et al.ll2005l . 120071 ) . In population simulations, the presence of a circumbinary disk 



can be at least as important as magnetic braking in the secular evolution of a CV, and can 
lead to completely dissolving the secondary star in a finite time, at an age ^ 3 times smaller 
than in an evolutionary model that does not include a circumbinary disk (jSpruit fc Taam 



20011 ). In our previous work, we have used Sp itzer mid-IR observations of several short 
orbital period (Pn rh = 1.3-1.8 hr) magnetic CVs (iHowell et al.ll2006l : iBrinkworth et al.ll2007l : 



Hoard et al.ll2007l ) to demonstrate an infrared excess caused by the presence of cool (T < 1000 



K), likely circumbinary, dust. In all cases, the total mass of dust required to reproduce the 
observed mid-IR spectral energy distributions (SEDs) was much smaller than the amount 
predicted to be necessary to have a significant effect on CV evolution. We now present 
the first analysis of Spitzer observations of a bright novalike (non-magnetic) CV at a longer 
orbital period, which offers a more likely environment for a substantial dust presence. 



2. Observations and Data Processing 



We observed V592 Cas with each of the three instruments on-board Spitzer (jWerner et al 



20041 ) during the G02 cycle. Details of the observations and data processing are given below. 



2.1. Infrared Array Camera 



We obtained photometric measurements in each of the four imaging channel^ of the 
Infrared Array Camera (IRAC; iFazio et al.ll2004l ) using 5 medium scale, cycling dithers of 
the 12-second frame time, which are available as standard options of the IRAC imaging 
Astronomical Observing Request (AOR) template. The observation (AOR identification 
number 14408960) was executed on 18 Aug 2005 UT as part of IRAC campaign 24. The 
data were processed through the 5])zteer Science Center (SSC) S14.0.0 pipeline to yield basic 
calibrated data (BCD) images. We applied the IRAC array location dependence correction 
to these BCDs, as described in the IRAC Data Handboolo. 

We did not perform photometry on a mosaicked image created from the BCDs becaus e 
of the poor results we found when attempting this in previous work (IBrinkworth et al.ll2007l ). 



^Channels 1-4 are referred to by corresponding wavelengths of 3.6, 4.5, 5.8, and 8.0 /im; however, see 
gSHHand Table [I 



'See http : //ssc . spitzer . caltech . edu/irac/dh/ 
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Instead, we performed aperture photometry on the individual corrected BCD images using 
the IRAF§ task QPHOT. We utihzed a 3-pixel radius aperture (1 pixel ~ 1.22"), with a 3-7 
pixel background annulus, which is one of the configurations for which aperture corrections 
are available in the IRAC Data Handbook (factors of 1.124, 1.127, 1.143, and 1.234 for 
channels 1-4, respectively). 

The individual IRAC channel 1 BCD images display a column artifact caused by a bright, 
saturated field star near the edge of the field of view. This column passes through the edge 
of the V592 Cas stellar profile, near the first diffraction ring. We investigated correcting 
for this artifact by interpolating across adjacent pixels. However, it became clear that the 
likely effect of this artifact is to decrease the nominal value of the channel 1 photometric 
measurement by less than la, so we did not include the correction to the final photometric 
measurement used in our SED of V592 Cas. The BCDs from the other IRAC channels do 
not contain this artifact because the bright star is not saturated. 

As described in the IRAC Data Handbook, we applied the pixel phase dependence 
correction to the channel 1 (3.6 /im) photometry. However, we did not perform the color 
correction other than to utilize the isophotal effective channel wavelengths (see Table [T]) 
during our subsequent interpretation of the data, which accounts for all but a few per cent 
of the color correction. The remaining effect of the color correction is folded into our IRAC 
systematic uncertainty budget - see below. We then performed a 6a clipping rejectiojlon the 
set of 5 individual BCD flux densities, which resulted in the rejection of only 1 measurement, 
from the channel 3 (5.8 /xm) data. Finally, the individual BCD photometry values for each 
channel were averaged together. 



2.1.1. Uncertainties in the IRAC Photometry 

For IRAC photometry, the total uncertainty budget includes several systematic terms: 
3% for the absolute gain calibration, 1% for repeatability, 3% for the absolute flux calibration 
of the IRAC calibration stars, and 1% for the color correction. These values are from 
version 3.0 of the IRAC Data Handbook. In addition, there is a scatter term, which we 



have determined as follows (based on the method described in iBrinkworth et al.ll2007l ). We 



performed additional aperture photometry on the BCDs, as described above, for all point 



■^The Image Reduction and Analysis Facility is maintained and distributed by the National Optical 
Astronomy Observatory. 

"^That is, a given flux density is compared to the average and standard deviation (cr) of all of the other 
flux density values. The measurement is rejected if it is more than 5cr away from the average. 
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sources (except V592 Cas) in the mutually overlapping regions of each IRAC channel in each 
of the 5 dithered exposures from both the on- and ofF-target fields of the two IRAC channel 
pairs (i.e., resulting in a total of 10 sets of photometry for each channel, comprised of up to 5 
measurements of each object in two different fields). We applied the 5a clipping rejection to 
the photometry, and kept all objects for which 3 or more (out of a maximum of 5 possible) 
measurements were available. This left a total of 543, 446, 112, and 62 objects in channels 
1-4, respectively. We then plotted, for each channel, the standard deviation of the individual 
flux densities for each object as a function of its average flux density. This results in well- 
behaved correlations between the scatter and flux density in each channel, which we show 
in Figured] for future use. We then determined what we call the "field-scatter" uncertainty 
for V592 Cas by averaging the points on these plots within ±20% of the mean flux density 
in each channel. This results in field-scatter terms of 0.05 (1.2%), 0.06 (2.1%), 0.08 (3.9%), 
and 0.08 (6.4%) mJy in channels 1-4, respectively. The field-scatter term is then added 
to the total systematic uncertainty budget described above. This process is more robust 
than simply using the "self-scatter" term (i.e., the standard deviation of the individual BCD 
photometry values for V592 Cas itself) since it is based on the results from a large number 
of objects and is not influenced by any intrinsic variability of the CV. 



2.2. Infrared Spectrograph Peak-Up Imaging Array 



We obtained a single photometric me asurement with the blue Peak-Up Imaging (PUI) 
array of the Infrared Spectrograph (IRS; iHouck et al.l |2004| ) using 30 small scale, cycling 
dithers of the 30-second ramp duration, which are available as standard options in the IRS 
PUI science mode AOR template. The observation (AOR identification number 14412032) 
was executed on 16 Jan 2006 UT as part of IRS campaign 28. The data were processed 
through the SSC S15.3.0 pipeline to yield basic calibrated data (BCD) images and a single 
post-BCD mosaic image. 

We performed aperture photometry on the post-BCD mosaic using the QPHOT task 
in IRAF. In order to avoid co ntamination from the relati vely bright nearby star 2MASS 
J00205373+5542192 (star 14 in iHenden fc Honevcuttlll997l . located 13" east), which is the 
only other object detected in the IRS-PUI mosaic, we utilized a 3-pixel radius aperture (1 
pixel = 1.8"), with the background computed from a mean over the entire image (excluding 
the two stars visible in the frame). The aperture correction for this configuration (factor of 
1.56) is available in the IRS Data Handboolj^. As with the IRAC data, we did not perform the 



http : //ssc . spitzer . caltech . edu/irs/dh/ 
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color correction, other than to utihze the isophotal channel wavelength (see Table [T]) during 
our subsequent interpretation of the data. The remaining effect of the color correction is 
folded into our IRS-PUI uncertainty budget (see below). 



2.2.1. Uncertainties in the IRS-PUI Photometry 

For IRS-PUI photometry, the total systematic uncertainty budget is the quadrature sum 
of 2% for the flux density zero point, 3% for repeatability, 3.4% for deviation from blackbody 
source SED in the aperture correction, and 5% for (reasonable) deviation from a fiat SED 
in the color correction. These values reflect upper limits from version 3.2 of the IRS Data 
Handbook. In addition, the total uncertainty includes a scatter term of 0.07 mJy (12.5%) 
obtained for the target aperture in the uncertainty ( "unc" ) image provided by the standard 
processing pipeline for the post-BCD mosaic. 

Because only two sources are visible in the small field-of-view of the IRS-PUI image, 
we could not repeat the field-scatter determination used with the IRAC data. We note that 
the average fiux density of V592 Cas obtained by performing photometry on each of the 30 
individual BCD images is 0.56 ± 0.12 mJy (where the uncertainty is the self-scatter term 
only), which is in agreement with the value measured from the post-BCD mosaic (see Table 
[T]). The larger scatter term in this case (which amounts to 21% of the measured fiux density) 
is likely a result of the fact that V592 Cas is very faint in any individual BCD, resulting in 
low S/N individual measurements. By comparison, the other source in the IRS-PUI field is 
brighter, with a mean fiux density of 1.92 mJy and a scatter uncertainty amounting to only 
7.5%. 



2.3. Multiband Imaging Photometer for Spitzer 



We obtained a single photometric measurement with the 24 fim array of the Multiband 
Imaging Photometer for Spitzer (MIPS; iRieke et al.l 120041 ) using 20 cycles of the 10-second 
exposure time small field photometry sequence, which are available as standard options in 
the MIPS pointed imaging AOR template. The observation (AOR identification number 
14412288) was executed on 21 Feb 2006 UT as part of MIPS campaign 29. The data were 
processed through the SSC S16.0.1 pipeline to yield basic calibrated data (BCD) images and 
a single post-BCD mosaic image. We re-mosaicked the data using MOPEX (version 16^ 



http : //ssc . spitzer . caltech . edu/postbcd/ onlinedocs /index . html 
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with overlap correction turned on, after first performing a self- calibration of the images; that 
is, we created a median image from all of the BCDs, normalized it to a value of 1.0, then 
divided all of the BCDs by that median image. This procedure removes persistent image 
artifacts that are not removed during the standard pipeline processing, such as long-term 
latents, low-level "jailbars" (that is, parallel linear artifacts produced following the read-out 
of a saturated source in the field), and gradients (see the MIPS Data HandboolJ^ for details 
of the self-calibration procedure and MIPS image artifacts). We resampled the pixel scale 
in this mosaic to 1.225 arcsec pixel"^ (i.e., oversampled by a factor of 2 compared to the 
standard post-BCD mosaic). 

We performed aperture photometry on our mosaic using the QPHOT task in IRAF. We 
utilized a 7-arcsecond radius aperture, with a 20-32-arcsecond background annulus. The 
aperture correction for this configuration (factor of 1.61) is available in the MIPS Data 
Handbook. We did not utihze the standard 7-13-arcsecond background annulus in order to 
exclude the star located 13" east of V592 Cas (see §2.2p . As with the data from the other 
instruments, we did not perform the color correction other than to utilize the isophotal 
channel wavelength (see Table [I]) during our subsequent interpretation of the data; the 
remaining effect of the color correction is folded into our MIPS uncertainty budget (see 
below). 



2.3.1. Uncertainties in the MIPS Photometry 

For MIPS photometry, the total systematic uncertainty budget is the quadrature sum 
of 4% for the absolute calibration, 0.4% for repeatability, and 3% for the color correction. 
These values reflect upper limits from version 3.3 of the MIPS Data Handbook. In addition, 
the total uncertainty includes a scatter term of 0.06 mJy (14%) obtained for the target 
aperture in the uncertainty ("munc") image provided by MOPEX for our mosaic. 

Because only a few bright sources are visible in the individual MIPS BCD images, we 
could not repeat the fleld-scatter determination used with the IRAC data, nor could we 
perform photometry for V592 Cas itself in the individual BCDs. We note that the average 
flux density of V592 Cas obtained by performing photometry on the standard post-BCD 
pipeline mosaic image is 0.45 ± 0.06 mJy (where the uncertainty is the scatter term only), 
which is consistent with the result from our improved mosaic (see Table [ p. The detailed 



charac terization of the properties of MIPS 24 /xm photometry presented in Ivon Braun et al. 



( 120081 ) suggests a somewhat lower fractional scatter uncertainty should be obtained for a 



^See http : //ssc . spitzer . caltech . edu/mips/ dh/ 



- 8 - 



source with the flux density of V592 Cas, in the range 7-13% (see Figure 6 in that work). 
However, the shghtly larger fractional scatter uncertainty that we observe in our MIPS 
photometry (14%) likely results from the fact that our single photometric measurement was 
obtained from a mosaic created from 280 individual BCD im ages (whereas each photometric 
measurement shown in Figure 6 of Ivon Braun et al.l |2008| was obtained from the average 
of 1224-2448 individual BCDs) coupled with the slightly higher 24 /zm background during 
the MIPS observations of V592 Cas (22.9 MJy sr-^) compared to GU Boo (19.7-21.0 MJy 
sr~^)§. With that in mind, our MIPS ph otometry results are con sistent with the expected 
level of photometric uncertainty found by Ivon Braun et al.l (120081 ) for a source of this mean 
brightness. 



3. Spectral Energy Distribution Model 
3.1. The Data 



Table [T] lists our observed Spitzer data for V592 Cas, as well as adopted flux densities 
at shorter wavelengths. The latter consists of two UV photometric measu rements from an 



Inter national Ultraviolet Explorer ( lUE) spectrum (as shown in Figure 1 of IWitherick et al. 



20031 : also see iGuinan fc Sionlll982l ). and ground-based V (see below) and 2MASS near-IR 
J, H, and Kg photometric measurements. 

When constructing an SED from these data, we do not have to be concerned about 
orbital variations since V592 Cas is a low inclination, non-eclipsing system. However, we do 
have to assume that non-coherent variability on the order of a couple tenths of a magnitude 
(i.e., normal CV "flickering") is present, as well as consider the possibility that there could 
be long-term trends in the brightness of V592 Cas. To aid in this assessment, in Table [2] we 
present a compilation of published optical and near-IR photometry of V592 Cas that is, to 
the best of our knowledge, complete!. 

The AAVSO measurement was obtained contemporaneously with our Spitzer observa- 
tions, and shows that V592 Cas was relatively cons tant in brightness during that time. The 
most recent l^-band measurement {{V) = 12.55; iKafka et al.l l2008l ). which was obtained 
^ 2 months prior to the IRAC observation, is somewhat brighter than historic observations 



^We note that one-half of the GU Boo MIPS data presented in von Braun et al. ( 20081 ) were obtained 
during the same MIPS instrument campaign, and on the same day, as our MIPS observation of V592 Cas. 

^Some photometry results, particularly the UBV measurements from Haue ( 1970l) . are repeated in mul- 
tiple references. In such cases, we have only presented the original reference source here. 
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of V592 Cas (V 



Haue 


1970; 


Huber et al. 


1998: 


Tavlor et al. 


1998) 



Zwitter fc Munaril (jl994j ) determined a color index of {V — R) = +0.11 for V592 Cas from 
their spectroscopic observations, so the AAVSO measurement corresponds to V ^ 12.4, 
which is consistent within the r a n^e of uncertainty and flickering variability of the l^-band 



measu rement from iKafka et al.l (120081 ). Consequently, we have adopted the iKafka et al. 
( 120081 ) value of = 12.55 for our SED. The only available near-IR photometry other than 
that from 2MASS appears to have suffered from calibration problems (see note in Table [2]), so 
we have used the 2MASS photometry values (in any case, these measurements do not appear 
to be significantly discrepant in the SED compared to the data at other wavelengths). The 
optical photometry was converted from magnitudes to flux densities using the zero points 
given in ICoxl ( 20001): for th e 2MA SS All Sky Point Source Catalog photometry we used the 



zero points in 



Cohen et all (120031 ). 



The flux density uncertainties for the 2MASS data are the la photometric uncertain- 
ties propagated through the corresponding magnitude to flux density conversions. For the 
ground-based opti cal data, the uncertainty reflects the observed range of V magnitude from 
Kafka et al.l (l2008l ) propagated through the flux density conversion. The flux density uncer- 
tainty on the UV measurements, which were read from a plot of the spectrum, is estimated 
at 10%. Derivation of the flux density uncertainties for the Spitzer data from each instru- 
ment is described in §2.1.H §2.2. H and §2.3.11 In all cases, the "uncertainties" on the filter 
wavelengths represent the widths of the photometric bands. 



3.1.1. Reddening and Distance to V592 Cas 



Nandy et al.l (119751 ) - among others - derive a correlation of the equivalent width of 



the strong UV interstellar absorption band at 2200 A with optical color excess, E[B — V) 
and the band profile shape with the nature and size of the grains producing the observed 
extinction. In this manner, iGuinan &: SionI (jl982l ) used the lUE spectru m of V592 Cas frora 
5 Dec 1981 UT to obtain a color excess (reddening) of E{B — V)= 0.3. iTaylor et al.l (Il998l ) 
used observations of the stationary (i.e., interstellar) Na D absorption in V592 Cas to derive 
a lower reddeni r ig of E(B — V) = 0.07-0.20. They adopt a "compromise" value with the 
Guinan fc SIotI Jl982h result of E{B - V) = 0.25. 



In the optical spectra of V592 Cas presented in iKafka et al.l (120081 ) ^ the Na D li n es are 
contaminated by Tucson city light emission, so we could not repeat the iTaylor et al.l (119981 ) 
measurement. However, we reanalyzed the archival spectrum of V592 C as, by applying 



the IR AF task deredden, which utilizes the empirical extinction law of ICardelli et al. 



(119891 ). Removal of the 2200 A absorption feature required a reddening of E{B — V)^ 0.22. 



- 10 - 



This is consistent with the iTaylor et al.l (119981 ) value. 



We used our reddening value with the diffuse interstellar extinction law from lFitzpatrick fc Massa 



(120071 ) to deredden the observed ffux densities of V592 Cas. B ecause there is some uncer - 
tainty associated with the extinction curve (e.g., see Figure 9 in iFitzpatrick fc Massal 120071 ). 
the dereddening process increases the uncertainties of the dereddened flux densities relative 
to the corresponding observed flux densities. The reddening correction is largest at short 
wavelengths and becomes increasingly negligible in the IR. Consequently, to avoid unneces- 
sarily degrading the accuracy of our observed flux densities, we did not apply the dereddening 
correction to any photometric measurement for which the fractional change produced by the 
correction would be smaller than 25% of the fractional standard deviation of the observed 
measurement {cri^/fu)- In practice, this means that the dereddening correction was not ap- 
phed to any of the Spitzer data - see Table [TJ Figure [2] (top) shows the dereddened SED for 
V592 Cas and the corresponding dereddened flux densities are listed in Table [H 



Harrison et al.l (l2004j ) have updated the My-Porb relation for dwarf novae in outburst 



(jWarneiill987l ) by using distances from Hubble Space Telescope Fine Guidance Sensor trigono- 
metric parallax measurements. They have also shown that this relation appears to be ap- 
plicable to novalike CVs as well. We used this relation, with extinction appropriate for our 
reddening value {Ay = 0.66 mag), to determine an expected (dereddened) absolute mag- 
nitude for V592 Cas of Mv,corr = 4.9 or My^nhs = 4.1, wh ere the former is the inclination- 
normalized value (see description in iHarrison et al.l |2004J and references therein) and the 
latter is the value that would be observed at the nominal inclinat ion of V592 Cas (i = 28°). 
Using V = 12.55, this gives a distance to V592 Cas of 360 pc. iTaylor et al.l (119981 ) made 
a similar calculation using the old calibration of the My-Porb relation, and found distances 
of 350 pc (if interstellar extinction is ignored), 240 pc (including extinction), or their final 
adopted distance of 330 pc (assuming that V592 Cas is "a bit more luminous than dwarf 
novae"). It is not clear if they made the correction for inclination, which would have been of 
the same amplitude and direction as that produced by their assumption that V592 Cas was 
more luminous than an outbursting dwarf nova. In addition, the revision to the calibration 
of the My-Porb relation would also have acted to make their distance larger by about 3%. 
We have scaled all of our SED model components (see §3.21) to a distance of 360 pc. 



3.2. The Model 



Our custom-built I R SED modeling code for CVs is described in lBrinkworth et al.l (120071 ) 
and iHoard et al.l (120071 ). We have continued to improve the code by implementing a model 
component to represent the accretion disk in non-magnetic CVs (described below). For this 
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work, we have used a model consisting of a WD (denoted with subscript "wd"), secondary 
star (subscript "ss"), optically thick steady state accretion disk (subscript "acd"), and opti- 
cally thin circumbinary dust disk (subscript "cbd") to compare to the observed SED of V592 
Cas. The system and model parameters are listed in Table [31 and Figure [2] (bottom) shows 
the SED data for V592 Cas with the nominal model. Figure [3] shows a to-scale diagram of 
V592 Cas based on the nominal model. We discuss some details of the model components 
and parameter choices below. 



3.2.1. The White Dwarf 



The WD is represented by a blackbody curve with T„d = 45, OOP K, which was cho- 



sen as a "typical " WD temperature for novalike CVs (e.g., lAraujo-Betancor et al.l 12003 



Hoard et al.ll2004l ). We have kept the WD temperature fixed at this value throughout our 
model calculations; for example, ignoring the possible role of the rate of mass accretion in 
heating the WD. In any case, as described below, even at this temperature the WD is only 
a minor contributor to the total system flux distribution longward of the ultraviolet (UV). 
If the temperature of th e WD in y 592 Cas is significantly different than this, then it is 
most likely to be smaller (jSiorull999l ). which would make the WD even less significant in the 
determination of a model that can reproduce the long wavelength observations. 

The WD is assumed to have the mean WD mass in CVs with Pnrh < 6 hr, Mwh = 0.75 Mr.) 



(iKniggd 120061). a nd a corresponding radius of -Rwd = 0.0106 i?© (IHamada fc Salpeteiill961 



Panel et al.ll2000l ). The projected surface area of the WD has been reduced by 7%, consistent 
with obscuration by an accretion dis k (see below) that extends outward from -Racd,in = 1-Rwd 
at the system inclination oi i = 28° (IHuber et al.lll998l ). 



3.2.2. The Secondary Star 



The empirical CV secondary star sequences of ISmith fc DhillonI (119981 ) and iKnigge 



(120061 ) predict a secondary star spectral type of M4.5±l V (their Equation 4) and M4±2 V 
(his Figure 9), respectively, for the 2.76 hr orbital period of V592 Cas. That is, regardless of 
the true nature of the secondary stars in CVs, which likely differ in internal structure from 
corresponding single main sequence stars (e.g., by being oversized and /o r slightly evo lved 
as a result of their history of mass loss and irradiation - see discussion in iKniggd |2006| and 
references therein) and/or by being non- uniformly heated by irradiation from the WD and 
accretion disk, the secondary star in V592 Cas has the spectroscopic appearance of a main 
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sequence star with spectral type in the range indicated above. As it turns out (see Figure [2]), 
the nominal secondary star contributes < 2.5% of the total flux density at all wavelengths, 
so the choice of secondary star template is robust against significant changes in spectral 
type or overall stellar temperature/luminosity that might be produced by evolutionary or 
irradiation effects (e.g., changing to an M2 template increases the maximum secondary star 
contribution to only ~ 4% of the total SED). 

Thus, to represent the secondary star, we have used an e r npirical template o f an M5.0 
V star (GJ1093) constructe d from optical jMermillio"dl llOsd Monet et al.l[2003h . 2MASS 
near- IR, and Spitzer IRAC dPatten et allbooeh photometry. The secondary star template 
was extended to wavelengths longer than 8 /xm by extrapolating the shorter wavelength data 
along a Rayleigh- Jeans tail (i.e., fi, oc A^^). The secondary st ar template was scaled fr om 
the trigonometric parallax distance of d = 7.76 pc for GJ1093 (lOppenheimer et al.ll200ll ) to 
the adopted distance of V592 Gas, d = 360 pc (see §3.1.11) . Avera ge mass, in itial radiu^, 
and temperature values for a star of this spectral type (compiled in iGoxl 120001 ) were used in 
the model. 



3.2.3. The Accretion Disk 



The accretion disk component assumes an optically thick steady state disk composed 
of concentric rings ern i tting as blackbodies, following the "standard model" prescription in 
Frank. King, fc Raind (120021 ). The accretion disk component is parameterized by inner and 
outer radii (-Racd,in and -Racd,out), height at the inner radius (/^acd), and a mass transfer rate 
from the secondary star (M). The inclination of the disk is assumed to be the same as the 
system inclination; higher inclination systems have smaller projected surface area, resulting 
in an overall decrease in the effective brightness of the disk. In the case of V592 Gas, the 
inclination is quite low {i = 28°) so obscuration of the WD by the disk (and vice- versa), 
although included in our model calculation, has an insignificant effect on the resultant SED. 

The general shape of the model accretion disk SED is similar to a slightly flattened 
and stretched blackbody function; that is, it displays a relatively broad peak with a rapid 
drop in brightness on the short wavelength side an d a gradual (Rayleigh-Jeans -like) decline 
on the long wavelength side (e.g., see Figure 20 in lFrank. King, fc Raind l2002l ). In general, 
increasing either the disk area (i.e., making -Racd.out larger) or the mass transfer rate causes 
the overall accretion disk SED to become brighter. Increasing -Racd.out also tends to shift the 



^"^The secondary star size and shape were recalculated within the Roche geometry of the binary in our 
model - see Figure El 
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peak of the accretion disk SED toward longer wavelengths, as more surface area is present 
in the cooler outer regions of the disk. Increasing M tends to shift the peak of the accretion 
disk SED toward shorter wavelengths, as the temperature in the inner disk increases. 

Our accretion disk models for V592 Gas were composed of 100 concentric rings, each 
of which contains 1441 azimuthal sections (corresponding to 0.25° per section). In addition, 
the outer edge of the accretion disk also contributes to the total disk SED. A wavelength- 
and temperature-dependent correction for limb darkening of the accretion disk is applied 
to each azim uthal section using li near limb darkening coefficients interpolated over the grid 



presented by IVan Hammd (ll993l l l^H In the case of a flat disk, this is a relatively simple 
procedure, as the disk face everywhere presents the same viewing angle (i.e., angle between 
the line-of-sight and a normal to the disk face), so only the disk edge is limb-darkened 
across a range of viewing angles. However, in the case of a flared disk (see below), which 
presents different viewing angles for different parts of the disk, limb darkening is calculated 
individually for each azimuthal section in each ring. The primary effect of limb darkening is 
to make the short wavelength end of the SED fainter relative to the long wavelength end. 



By default, the standard model radial temperature profile from iFrank. King, fc Raine 



( 120021 ) (their Equation 5.41) is used to determine the temperature in each disk ring. It is 
possible for the standard model radial temperature profile (T oc r~"' with 7 = 0.75) to be 
modified by irradiation; that is, by reabsorption in the outer disk of fiux from the central 
object and/or inner disk. Irradiation is primarily important in fiared disks (i.e., disks whose 
height increases with radius), which can more easily intercept luminosity from the WD and 
inner disk. In order to acc ount for this po s sibilit y, we have provided an option to utilize the 
irradiation prescription of lOrosz fc Wadd (120031 ). in which the radial temperature profile of 
the outer accretion disk (at i?acd > -Rcrit, where Rent = ((7 — 87out)/(6 — 87out))^ -Rwd) is 
allowed to have a shallower temperature gradient than in the standard model (7out < 0.75). 
At the same time, the height of the d isk for radii larger than Rent is assumed to increase as 
R^/^ (equat ion 2.51b in IWarnerlll995l ). The outer edge of the disk, which does not "see" the 
WD and inner disk, is always assumed to have a non-irradiated temperature appropriate for 
the standard model radial temperature profile. 



^^We als o tested a wavelength-dependent quadratic limb darkening law based on observations of the Sun 
but this did not produce any significant difference compared to the linear limb darkening law. 
Limb darkening corrections were not applied to the stellar components. In the case of the secondary star, 
this is because the component SED is based on empirical templates which are intrinsically limb -darkened. In 



the c ase of gray atmosphere limb darkening, which may be most ap propriate for the W D (e.g.. lWarner et al 



19711 ). the limb darkening correction is equivalent to a small (< 5%; 



Davis et al.ll2000l ) change in the assumed 



radius of the WD, which is already uncertain by at least this amount. In addition, limb darkening is not 
applied to the circumbinary dust disk, since it is treated as optically thin. 
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As found in lOrosz &: Wadd (|20CI3l) an d other studies of the effect of irradiation on 
model accretion disk SEDs (e.g., Wadel 1988 ). we found that irradiation (hke hmb darkening) 



primarily affects the short wavelength end of the SED. However, irradiation (which effectively 
raises the temperature at a given radius over the non-irradiated case) tends to make the short 
wavelength end of the SED brighter relative to the long wavelength end, which somewhat 
counteracts the effect of limb darkening. Changing (increasing) the mass transfer rate can 
often compensate for the combined effects of limb darkening and irradiation in the model 
accretion disk SED compared to the fiat disk (i.e., no irradiation) case. However, for small 
values of 7out, the model SED becomes too faint at the short wavelength end, such that 
increasing M enough to supply the "missing" UV flux causes the accretion disk to exceed 
the observed SED at longer (e.g., IRAC bandpass) wavelengths. For V592 Cas, this occurs 
for 7out ^ 0.60, which requires M > 1.8 x 1O~^M0 yr~^ to reproduce the short wavelength 
end of the SED, but then exceeds (with the accretion disk component alone) the observed 
flux densities in the 2MASS and IRAC bands. Figure [2] shows SED models for V592 Cas 
using either a limb-darkened fiat accretion disk or a limb-darkened flared accretion disk with 
irradiation (7out = 0.7). There is no significant difference in the resultant total model SEDs. 
The corresponding model parameters for each case are given in Table [31 In the remaining 
analysis, we utilize the flared accretion disk model. 

It is not surprising that the accretion disk of a nearly face-on novalike CV is very bright. 
Indeed, the nominal accretion disk component for V592 Cas dominates at all wavelengths 
from the optical through near-IR, and even into the mid-IR (IRAC) range, accounting for 
all but < 10% of the total flux density in the wavelength range 0.5-5 /ino. This makes 
accounting for the flux density that is in excess of the accretion disk model component at 
long wavelengths (A > 8 /im) the principle challenge to reproducing the observed SED of 



^^Many CVs are inferred to have a bright spot on the edge of the accretion disk at the site of the 
initial accretion stream impact. In principle, this could introduce another component into the model. We 
experimented with representing the bright spot in our V592 Cas SED model by using an additional blackbody 
component with temperature of 20,000 K (i.e., peaking in the near-UV) and projected emitting area up to 
several times the projected area of the WD. However, the result is a contribution of only < 1% of the total 
flux density at all wavelengths, and < 0.5% at all wavelengths longer than 1 fjja. Since the flux density 
uncertainties at the blue end of the SED range from 10-20% from the optical to UV (mainly due to the 
dereddening correction), the contribution of a disk bright spot is below the level of detect ability. In a higher 
inclination CV, in which the relative contribution from the acccretion disk face is reduced, a bright spot 
component could be an important addition to the SED model. Additional components, such as a boundary 
layer, optically thin disk corona or wind, etc., could also contribute to the accretion disk SED, but are 
not considered here because they would primarily affect only the UV end of the spectrum, have relatively 
sparse empirical and theoretical bases in the mid-IR, and would exceed the ability of the current data set to 
constrain. 
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V592 Cas. 

Figure H] shows the effect of changing -Racd.out (top panel) or M (bottom panel) in the 
nominal flared accretion disk component. It is clear from this figure that adjusting either 
parameter individually can not account for the flux density excess at long wavelengths. In 
the case of -Racd.out, the maximum allowed value is ~ 45i?wd, corresponding to the size of the 
WD Roche lobe, and even this is insufficient to reproduce the 16- and 24-/im measurements. 
Increasing M also fails - when the model is bright enough to match the 16-/im measurement, 
then it is not bright enough to match at 24 /im, and when it would be bright enough to 
match the 24-/im measurement, then it would be too bright at 16 /xm. The latter case 
would require an exceptionally high mass transfer rate, in excess of 10~® yr~^, which is 
not plausible. Even worse, in order to match even the 16-/im measurement by increasing 
either parameter (even disregarding the firm upper limit to -Racd.out), the model flux density 
greatly exceeds the observed values at shorter wavelengths. Changing both parameters in 
unison (e.g., increasing -Racd,out while decreasing M) also does not produce a viable solution, 
since this does not change the slope of the long wavelength end of the accretion disk model 
component - it will never match both the 16- and 24-/im measurements at the same time, 
and will always be too bright at shorter wavelengths when it is forced to match either of 
the longer wavelength measurements. We are left with the conclusion that, in the context of 
the standard steady state accretion disk model (with or without irradiation), an additional 
component must be present to explain the observed 16- and 24-/im data - this is discussed 
in §2231 

The nominal accretion disk outer radius is -Racd,out = 35i?wd, which is equivalent to 
about 75% of the WD Roche lobe radius. This is consistent with estimates for the sizes of 
accretion disks in novalike CVs, which, at the mass ratio of V592 Cas, range from ~ 0.3-Rli 
from the theory for a ze ro viscosity disk to ~ 0.4-0 . 8-Rt,i from disk reconstru ctions based 



on eclipse observations (IRitterl llQSOl : ISulkanen et al.l Il98ll : iRutten et al.lll992l ). The mass 
transfer rate, M ^ 1 x 10~^ Mq yr~^ is in the middle of the ^ 10~^-10~^ M p, yr~^ range that 



is often cited as characteristic of novalike CVs (IWarnerlll995l ). iTaylor et al.l (Il998l ) estimated 



a mass transfer rate of 9 x 10~^ Mq yr~^ for V592 Cas from a power law fit to an optical- 
near-IR SED, which is quite close to our result (we note that if we utilize a fiat standard 
model accretion disk with no irradiation and no limb darkening, the n our best match t o the 
observed SED is obtained with M = 9 x 10^^ Mq yr~^, as found by lTaylor et al.lll998l ). 



The orbital period of V592 Cas places it at the upper end, or well inside, the CV 
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extrapolating between the expect ed mass transfer rates at either end of the period gap (e.g., 



Warner! I1995I : iHowell et aLll200ll ). we would expect that at the orbital period of V592 Cas, 
in the absence of whatever mechanism causes the gap, i t should have M ~ 1O~^M0 yr~^, 
an order of magnitude smaller than inferred here and in iTaylor et al.l (119981 ). However, the 

iod end of the gap ranges from ~ 10~^''-10~* Mq yr~^ 
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the period gap but still accreting at a high rate implies that it is not a run of the mill CV. It 
has either recently evolved into the gap and is in the process of shutting down mass transfer, 
or is a relatively young CV that only recently initiated mass transfer at close to its current 
orbital periocO. 



3.2.4- The Circumhinary Dust Disk 



In order to supply the requisite additional flux density at the longest mid-IR wave- 
lengths, we utilize a circumbinary dust disk component. This model component is calcu- 
lated as described in iHoard et al.l (120071 ). under the assumption that the disk is optically 
thin, and composed of spherical grains with radius of 1 /xm that (re)radiate as blackbodies. 
The problem of determining the inner boundary condition of temperature in the presence of 
m ultiple irradiating b odies (i.e., the WD, secondary star, and accretion disk), as described 



m 



Brinkworth et al.l (120071 ). is also an issue for V592 Cas. Because the WD in V592 Cas 



is relatively hot, and both the secondary star and accretion disk are large, we expect that 
the maximum radius at which the dust temperature would be above 1000 K (i.e., in the 
1000-2000 K range for dust sublimation) is on the order of 200-300 R^^ in V592 Cafl 
This is well beyond the tidal truncation limit of ~ 1.5 times the binary separation (about 
157i?wd for V592 Cas), which we have assumed for the inner edge radius of the circumbinary 
di sk in previous works. Con sequently, we have used the "equivalent star" prescription given 



m 



Brinkworth et al.l (120071 ) to set the temperature at the inner edge of the circumbinary 



disk for a given radius. We note that essentially identical circumbinary disk SEDs can be 
obtained using different combinations of inner and outer radii by adjusting the total mass 
in the disk (e.g., we can obtain an SED identical to the one shown in Figure [2] by forcing 



^■^A third possibility is that V592 Cas is an old CV which is evolving to long er orbital period out of the gap 
due to the evolutionary effect of a massive circumbinary disk (as discussed in lWillems et al.ll2007l) : however, 
in light of the low dust mass that we have found - see ij3.2.4l and i )4.4H4.5l - we consider this scenario unlikely. 
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This also likely rules out a substantial circumstellar dust ring around the WD - as discussed for the 



short orbital period dwarf nova WZ Sge in I Howell et al 
excess in V592 Cas. 



(200. 



3) - as the source of the long wavelength IR 
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the circumbinary disk inner edge radius to be at the tidal truncation radius, if the total 
mass of dust in the disk is increased by a factor of ~ 5). Similarly, the effect of different 
prescriptions for the radial temperature profile for the dust can be offset by adjusting the 
inner /outer radii and/or total mass of the circumbinary disk. In the context of reproducing 
the observed data, the important parameter is the temperature at the inner edge of the disk 
- it cannot be significantly higher than ~ 500 K, or the circumbinary disk component SED 
becomes too bright at short wavelengths to match the observations. 



4. Discussion 



4.1. Excluding Bremsstrahlung Emission 



Howell et al.l (120061 ) discuss several general reasons why bremsstrahlung, or free- free, 
emission in an outflowing w ind can be excluded as a significant source of mid-IR emission 
in CVs. iDubus et al.l (120041 ) present quantitative observational tests that would confirm or 
refute IR emission in CVs as being caused by bremsstrahlung emission. The first of these 
relates to the expected fi, oc u^'^ spectrum of bremsstrahlung emission. At 16 /im, the IR ex- 
cess in V592 Cas (i.e., the remaining flux after the wd-|-ss+acd model components have been 
subtracted) is 0.33 mJy; at 24 fim, the excess is 0.30 mJy. This is a ratio of fu,24:/ fu,i6 = 0.91, 
whereas if the IR excess followed a spectral law consistent with bremsstrahlung emission, 
then we should expect the flux density at 24 /xm to be ~ 0.78 x the flux density at 16 fim. 

The second test relates to the wind mass loss rate required to e xplain the observed 
level of IR excess via bremsstra hlung emission. Using E quation 4 from iDubus et al.l (120041 ) 
(which is bas ed on the work of IWright fc Barlowl 119751 ) with an assumed wind velocity of 
2000 km s~^ (IKafka et al.l l2008l found wind velocities reaching up to 5000 km s~^ in V592 
Cas), the observed IR excesses at 16 and 24 fim would require wind mass loss rates of ~ 4.0- 
5.2 X 10"^ M0 yr^^. This is equivalent to ~25-35% of the entire inferred mass transfer rate 
from the secondary star in V592 Cas and, as such, is improbably high for the wind mass loss 
rate. 



4.2. No Dust Emission at Short Wavelengths 



Belle et al.l ( 120041 ) concluded that there is no evidence for the presence of a circumbinary 



dust disk in the IR SED for V592 Cas; however, their data spanned only the JHKL' bands 
(i.e., ~ 1.2-3.8 yum). From our data, it is clear that the IR excess in V592 Cas is not 
readily apparent until well past 5 /im. For example, at wavelengths < 5 /im, the wd+ss+acd 
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components account for > 90% of the observed flux density, wliereas at wavelengths > 8 
/im, these components account for < 70% of the observed flux density. At L\ only 3% of 
the observed flux density is not accounted for by the wd+ss+acd components. In addition, 



Belle et al.l ( 12004 ) were looking for a much hotter circumbinary disk component, with T 



3000 K, whereas our observations show that the maximum temperature of the dust is T ^ 500 
K. We want to emphasize that this is not just a case in which a bright, near-face-on accretion 
disk is "masking" fainter dust emission at short wavelengths. Dust at hotter temperatures 
in V592 Cas would have been detectable at A < 5 /xm, even against the comparatively bright 
contribution of the accretion disk at those wavelengths, if the long wavelength end of the 
dust distribution was bright enough to reproduce the observed 16- and 24-/im measurements. 



4.3. Dissipation and Formation of the Dust 



One question that has plagued a full understanding of the presence of circumbinary 
disks in CVs is the origin of the dust. Possible scenarios include: a remnant of the common 
envelope phase of the CVs early evolution, formation through tidal disruption of minor 
planets or comets that survived from the progenitor solar system (as assumed for the origin 
of circumstellar dust disks around single WDs), creation via mass outflows in winds and nova 
outbursts from the inner binary (the m echanism assumed fo r the formation of circumbinary 
dust disks in the theoretical models of IWillems et al.l 120071 ). or some combination of these 
processes. Regardless of the exact origin of the dust, since the calculation of the structure of 
the circumbinary disk is predicated on the assumption that there is a region in the CV system 
in which temperatures are too high for dust to exist, it seems likely that the circumbinary 
disks cannot be entirely static: dust will be sublimated at some rate at the inner edge as 
grains stray too close to the inner binary, and another trickle of dust will be lost at the outer 
edge of the circumbinary disk as grains stray too far away from the gravitational influence 
of the inner binary. So, it might be more appropriate to wonder about not only the origin 
of the dust, but its longevity as well. In V592 Cas, at least, there is an obvious mechanism 
to provide a continu ing supply of raw mate ri als for the circumbi nary disk, namely, the w ind 
outflow discussed in Iwitherick et al.l J2003h . IPrinia et al.l J2004h . and lKafka eTaP J2008h . 



The total circumbinary dust disk mass in V592 Cas amounts to only ~ 10~^^ M©, which 
is a small fraction (~ 0.01%) of the total annual mass transfer budget in this system. If the 
total mass flux present in the wind outflow from V592 Cas is a tenth of a percent of the total 
mass transfer budget, and even if the efficiency of capturing material from the wind into 
the circumbinary disk is a similarly small fraction, then this mechanism can still completely 
replenish the circumbinary disk on a timescale of ~ 100 yr. Another way to look at this 
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is that as long as the mass loss rate from the circumbinary disk (e.g., dissipation through 
sublimation and/or escape from the system) is smaller than the rate of dust formation from 
the wind, then the wind can maintain the circumbinary disk at a constant mass. Using 
the 0.1% efficiencies for the wind outflow and dust formation rates, as suggested above, 
the circumbinary disk mass will be maintained as long as < 1% of the total circumbinary 
disk mass is lost per year. This result also implies that in V592 Cas a much more massive 
circumbinary disk probably could not be maintained by the wind outflow alone. 



4.4. Total Mass of Circumbinary Dust 



The total mass in our model circumbinary dust disk is Mchd ~ 10^^ g. Reasonable 
variations in the inner radius of the circumbinary disk that would be required to yield 
the necessary boundary condition of Tcbd,in ~ 500 K in the presence of different irradiation 
prescriptions (i.e., treatment of multiple heating sources in the inner binary) can be balanced 
against changing the total mass in the circumbinary dust disk to produce an SED that is 
indistinguishable from the one presented here. In addition, the photometric uncertainties 
also provide some leeway for changing the dust mass while still adequately reproducing the 
observed SED. However, these variations in mass are relatively small (i.e., factor of < 5). 
The mass of the circumbinary disk in V592 Cas exceeds those found in several of our past 
i nvestigations of short o rbital period magnetic CVs, which were on the order of 10^^-10^^ g 
(IBrinkworth et al.l 120071 ). It is comparable to the high end of the range of m ass estimates 



for t he circumbinary disk in the short orbital period, magnetic CV EF Eri ( iHoard et al. 



20071 ). However, as found in all of these previous studies, the total mass of dust is many 
orders of magnitude smaller than the 10^^-10^^ g predicted to be necessary for circumbinary 
disks to serve as a ngular momenturn loss mechanisrns tha t can signiflcantly a ffect the secular 



evolut ion of CVs (ITaam et al.ll2003l : IWillems et al.l 120071 ). The prediction by IWillems et al. 



(120071 ) that the circumbinary dust disks in (disk-accreting) CVs should only be detectable 
at wavelengths longer than 10 fim is borne out by our observations for the case of V592 Cas. 

The outer radius of the circumbinary disk has been somewhat arbitrarily set at 15,000 
i?„d, giving an outer temperature of 50 K, comparable to what one might expect as the disk 
merges with the local ISM. The model SED at the observed wavelengths is not sensitive to 
decreasing this outer radius by more than a factor of 2, corresponding to an outer temperature 
of ~ 100 K. However, in order to continue to match the observed mid-IR photometry, 
the total mass of dust in the disk must be decreased as the outer radius decreases - this 
exacerbates the low mass discrepancy of the circumbinary disks compared to the expectation 
from CV evolution theory. Unfortunately, applying the reverse of this scenario also cannot 



- 20 - 



resolve the low mass problem: increasing the outer radius of the circumbinary disk to as much 
as 100,000 i?wd (corresponding to a rather low temperature of only 12 K) only increases the 
requisite total dust mass by a factor of ~ 7. 



4.5. Excluding an Optically Thick Circumbinary Disk 



The evolutionary calculations of iTaam et al.l (120031 ) and IWillems et al.l (120071 ) assume 
that the circumbinary disks are massive (> 10^^ g), hot (Tcbd,in ~ 1500-3000 K), and optically 
thick. If the circumbinary disks in CVs do not have these characteristics, then we should 
not expect them t o act as effective angula r mon aentum loss mechan isms during the course 
of CV evolution. iBrinkworth et al.l (120071 ) and iHoard et al.l (120071 ) discussed the reasons 
for preferring an optically thin circumbinary disk model over an optically thick one when 
attempting to match the observed IR excesses of magnetic CVs. For the sake of completeness, 
and because it would offer the simplest method for substantially increasing the dust mass 
in V592 Cas compared to the optically thin case, we have also explored the optically thick 
circumbinary disk case. We can produce a model circumbinary disk SED very similar to 
that shown in Figure [2] using an optically thick treatment. However, we must force the 
circumbinary disk to have an inner edge temperature of Tcbd.in ~ 500 K at the tidal truncation 
radius - this is a factor of ~ 3-4 cooler than the expected temperature at this radius due to 
irradiation from the WD. Starting the circumbinary disk at a larger radius where the expected 
temperature is closer to ~ 500 K does not work because the correspondingly larger surface 
area of the circumbinary disk makes the SED component too bright to match the observed 
mid-lR data. In addition, the circumbinary disk must be truncated at an outer radius 
corresponding to a temperature of ~ 250 K, much warmer than the ambient temperature 
of the ISM, or the model SED component is too bright at the long wavelength end of the 
observed range. These inconsistencies again suggest that the optically thick case is not 
appropriate for matching the shape of the observed IR excess in V592 Cas. Thus, we are led 
to conclude that the results from the optically thin case are more likely to be representative 
of the distribution of dust. 



4.6. An Initially Massive Circumbinary Disk? 

Considering the low dust masses that we have found in several CVs (now including 
a high M novalike CV with a wind; that is, a good candidate to have a large amount of 
circumbinary material), we wonder if an initially massive circumbinary disk (perhaps formed 
during the common envelope phase of the CVs early history) could evolve over time into 
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a low mass remnant that we see now. To begin to address this question, we have used a 
"toy" model to explore the general behavior of an initially massive circumbinary dust disk 
in the presence of dissipation at the disk edges and contribution from formation of new 
dust using material in the wind from the inner binary, as described in §4.31 This model is 
purely parametric, and does not consider the precise physics of how a particular rate of dust 
dissipation or formation might arise. The total mass of circumbinary dust as a function of 
time (years) is simply calculated under the assumption that a constant small fraction (ei) of 
the total dust mass is lost each year through dissipation. At the same time, a small fraction 
{62) of the mass transfer rate in the CV (M) is carried away in a wind, of which another 
small fraction {e^) forms dust and contributes to the mass of the circumbinary disk. That 
is, 

M,bd,i = (1 - £l)M,bd,i-l + ^2^3^, (1) 

and this equation is evaluated over i = 1-10® yr, with initial mass of at least 10^® g at 2 = 0. 
We assume that log(M) oc — log(i) (i.e., the mass transfer rate decreases linearly in log- log 
space), with an initial value of Minitiai = 10^ ''•^M© yr~^, and evolve the system over 10® yr 
to a final mass transfer rate of 10~® Mq yr^^ or smaller. This evolving mass transfer rate, 
to first order, m imics the expected decrease in M during the evolution of a CV above the 



period gap (e.g., IShafterlll992l : IWarnerlll995l : iHowell et al.ll200ll ) 



Several examples of the outcome of this calculation are shown in Figure O In gen- 
eral, the dust mass follows a similar evolution in all cases, over a wide range of values for 
-^finah ^cbd,i=0) ^1) ^2, a^d £3. It first goes through a period of nearly exponential decline 
during which the rate of mass contributed to the circumbinary dust as new dust forms is 
insignificantly small compared to the total mass of dust, and the losses from the circumbi- 
nary dust (parameterized by Si) dominate. At some critical time, Tcrit (marked by a sharp 
"elbow" in the plot of log(time) vs. log (dust mass)), the rate of mass being contributed to 
the circumbinary dust by the wind (parameterized by 62 and £3) becomes comparable to the 
rate of mass being lost, and the total dust mass becomes relatively stable for the remaining 
evolution of the system (showing at most only a gradual decline by less than an order of 
magnitude over many tens of millions of years). In general, the value of r^rit and the final 
dust mass both increase as ei is decreased because mass is lost from the dust at a lower rate 
and, consequently, it takes longer for the contribution from the wind to become significant 
in comparison to the rate of dust dissipation and the total mass of dust. The value of Xcrit 
also increases, to a lesser extent, when the initial dust mass is larger, because there is more 
dust to dissipate before the contribution from the wind becomes significant; however, if all 
other parameters are fixed, the evolution of the dust mass after Tcrit is identical regardless 
of the value of Mcbd,i=o- 

Because £2 and are used as a product in Equation [H it does not matter for the 
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sake of calculation which parameter is changed; however, they do have different physical 
interpretations, as defined above. Changing the product of €2 and £3 has only a small effect 
on Tcrit for a given Si. This mainly affects the final total mass of dust, with the mass 
increasing as the product £2^3 increases, which causes the rate of dust formation out of the 
CV wind to increase and, consequently, become significant earlier compared to the total dust 
mass and rate of dust dissipation. 

Panel (a) in Figure shows the result of varying the circumbinary dust dissipation rate 
efficiency, Si = 0.01, 0.001, 0.0001, with the wind-fueled dust formation rate efficiencies fixed 
at 62 = €3 = 0.001. The result of this is to cause Tcrit to vary from ~ 10^-10^ yr and the 
final dust mass to vary from ~ 2 x 10^°-2 x 10^^ g as ei decreases. By comparison with the 
other cases (discussed below), it is clear that the value of ei has the largest effect on Tcrit. 

Panel (b) in Figure O shows the result of varying the product of the wind-fueled dust 
formation rate efficiencies, 8283 = 0.0001,0.00001,0.000001, with the circumbinary dust 
dissipation rate efficiency fixed at ei = 0.001. In this scenario, the value of Tcrit is relatively 
unchanged at ~ 10^ yr, but the final dust mass varies from ~ 2 x 10^^-2 x 10^^ g as £2^3 
increases. In both panels (a) and (b) the initial dust mass and final mass transfer rate were 
fixed at Mcbd,i=o = 10^^ g and Mgnai = 1O~^M0 yr~\ respectively. 

Panel (c) in Figure [5] shows the result of varying the final mass transfer rate, Mgnai = 
10^^, 10~^, 1O~^°M0 yr~^, with the dust dissipation and formation efficiences fixed at ei = 
^2 = ^3 = 0.001. For the dust mass evolution up to Tcrit, the result is very similar to the 
case shown in panel (b). After Tcrit, the total dust mass declines more slowly for larger 
values of Mfinai, since there is a correspondingly larger wind outflow from which new dust 
can form. The case for Mfinai = 10~^ Mq yr~^ corresponds approximately to V592 Cas at 
the current time, and shows the least change in total dust mass in the post-Tcrit evolution of 
the CV. The case for Mfinai = 10~^ Mq yr~^ might be considered more "typical" of CVs at 
this orbital period, or could also represent the future of V592 Cas after continued evolution 
further reduces its mass transfer rate. 

Finally, panel (d) in Figure E] shows the result of changing the initial dust mass, 
Mcbd,i=o = 10^*^, 10^°, lO^'' g, with the dust dissipation and formation efficiences fixed at 
ei = 62 = 63 = 0.001 and the final mass transfer rate fixed at Mfinai = 1O~^M0 yr~^. We 
had initially chosen 10^*^ g of dust as a lower limit starting point of interest; obviously, the 
case in which the initial mass of dust is already smaller than predicted to be necessary to 
affect CV evolution is of little interest. However, if the initial dust originates, for example, 
during the common envelope phase, then there could potentially be more than this lower 
limit case. However, as panel (d) shows, even an initial mass as much as 10^ times larger 
than the lower limit results in identical post-Tcrit conditions and evolution, with only a small 



-23- 



increase in Tcrit as Mcbd,i=o increases. 

In general, we find that this simple toy model suggests that even an initially very 
massive (>> 10^^ g) circumbinary dust disk present at the birth of a CV will rapidly evolve 
(in only 10^-10^ yr) into a disk with much lower mass (~ lO^'^-lO^^ g). The dust disk mass 
is then approximately stable over the remaining evolution of the CV. The only requirements 
to produce this scenario are that the circumbinary dust experiences ongoing low rates of 
dissipation (e.g., via sublimation and escape from the system's gravity) and formation (e.g., 
from material carried outward in a wind from the inner binary). If a massive circumbinary 
dust disk is a common by-product of the common envelope phase of CV formation, then 
such a scenario might explain why all of the CVs so far observed with Spitzer appear to 
have only small amounts of circumbinary dust. We are encouraged by the relatively close 
agreement between our inferred mass of dust in V592 Cas and the toy model prediction for 
the final mass of circumbinary dust in a high mass transfer rate CV that has evolved to the 
long period end of the period gap. 

To first order, the initial rapid decrease in the mass of circumbinary dust would seem 
to imply that circumbinary dust disks do not play a role in the secular evolution of the CV 
population. However, this begs a new question: is the presence of a massive circumbinary 
disk for a short time early in the evolutionary history of a CV sufficient to have a lasting 
impact on the secular evolution of the CV? That is, could such a scenario result in an early 
episode of extra angular momentum loss of sufficient magnitude to alter the characteristics 
at the far end of the evolutionary sequence of the CV population compared to predictions 
utilizing angular momentum loss via magnetic braking and gravitational radiation alone? 
The answer to that question is beyond the scope of this paper. 

5. Conclusions 

We have found compelling observational evidence for the existence of IR excess in the 
near-facc-on, novalikc CV V592 Cas. This result stems primarily from our new Spitzer IRS- 
PUI and MIPS observations at 16 and 24 /xm, respectively. We have provided a representative 
model that reproduces the observed (dereddened) SED using a combination of white dwarf, 
secondary star, steady state accretion disk, and circumbinary dust disk components. At 
short wavelengths (A < 5 //m), the SED of V592 Cas is dominated by the steady state 
accretion disk and, as found in past studies of this system, there is no indication of the 
presence of circumbinary material until the SED at longer wavelengths is considered. The 
total mass in the circumbinary dust disk of V592 Cas, consistent with recent studies of IR 
excess in magnetic CVs, is apparently too small to have a significant effect on the secular 
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evolution of the CV. This is the first estimate for the mass of circumbinary dust in a non- 
magnetic, disk-accreting CV at an orbital period (~ 166 min) significantly longer than the 
observed period minimum 80 min), and demonstrates that the low dust mass results from 
our previously studied systems are not solely a result of the fact that they contain strongly 
magnetic WDs. 

The problem of accounting for the origin and/or longevity of the circumbinary dust 
is mitigated in the case of V592 Cas by the presence of a wind outflow, which can easily 
provide the necessary raw materials to replenish the circumbinary disk on short timescales 
or maintain it at a relatively constant mass. An implication of this result is that CVs 
without significant wind outflows might be expected to contain little or no circumbinary 
dust. However, our previous Spitzer observations of magnetic CVs show that dust can be 
present even in systems that lack a strong wind. So, there must be more than one route to 
the presence of circumbinary dust in CVs. 

We have made a simple calculation that shows that an initially massive circumbinary 
dust disk, perhaps formed during the common envelope phase of the binary's early evolution, 
will rapidly evolve to a relatively stable low mass in the presence of ongoing low rates of 
dissipation of dust in the existing circumbinary material and formation of new dust (for 

example, fueled by matter carried outward in a wind from the inner binary). It remains to 
be demonstrated in detail whether or not such evolving circumbinary disks can still play a 
role in the secular evolution of the CV population. 

This work is based in part on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology, under a 
contract with the National Aeronautics and Space Administration (NASA). Support for 
this work was provided by NASA. This work makes use of data products from the Two 
Micron All Sky Survey, which is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/Caltech, funded by NASA and the NSF. This 
research has made use of the SIMBAD database, operated at CDS, Strasbourg, Prance. 
We acknowledge with thanks the variable star observations from the AAVSO International 
Database contributed by observers worldwide and used in this research. 
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Table 1. Spectral energy distribution data 



Flux Density: 



Band 


Wavelength 


Observed 


Dereddened^ 


Date Obtained 






(microns) 


(mJy) 


(mJy) 




(UT) 


UV 





1450 


19.99 ± 2.00 


120 ±25 


05 


Dec 


1981'^ 


UV 





1800 


20.53 ± 2.05 


105 ± 19 


05 


Dec 


1981'" 


V 





55 ±0.045 


3fi 1 1 +-^-36 
"^"■^^-3.65 


66.6 ±7.0 


08 


Oct 


2005 


2MASS-J 


1 


r,or+0.125 

^'-''-'-0.115 


19.27 ±0.53 


22.28 ± 1.85 


28 


Sep 


2000 


2MASS-H 


1 




12.821°;!^ 


13.95 ± 1.32 


28 


Sep 


2000 


2MASS-K, 


2 




8.88 ±0.25 


9.35 ±0.91 


28 


Sep 


2000 


IRAC-1 


3 


550 ±0.375 


4.15 ±0.19 




18 


Aug 


2005 


IRAC-2 


4 


493 ± 0.508 


2.92 ±0.14 




18 


Aug 


2005 


IRAC-3 


5 


731 ±0.713 


1.98 ±0.12 




18 


Aug 


2005 


IRAC-4 


7 


872 ± 1.453 


1.22 ±0.10 




18 


Aug 


2005 


IRS-PUI-blue 




0.56 ±0.08 




16 


Jan 


2006 


MIPS-24 




0.43 ±0.07 




21 


Feb 


2006 



^Flux densities d e redde ned using the UV-optical-IR extinction law from 



Fitzpatrick fc Massal (|2007l) with E{B - V) = 0.22. See text for details. 



''Note that lWitherick et al.l (|2003r ) mistakenly give the date of the /J7i? spectrum 
from which these measurements were taken as 1984. 



Table 2. Photometry of V592 Cas 



JD 


Year 


U 


B 


V 


R J 


H 


Ks 


Reference 


c. 2438761 


c. 1965 


12.07: 


12.86: 


12.79: 








Haug (19701 


2449394.20-2449401.28 


1993.112-1994.131 




12.5-12.7 










Taylor et al. fl998) 


2450302.68-2450303.72 


1996.600-1996.603 






12.8(1) 








Huber et al. ( 1998) 


2450357.84 


1996.751 








12.8(2)=' 


12.2(2)=^ 


12.3(2)=' 


Huber et al. f 1998) 


2450741.64-2450836.62 


1997.803-1998.063 






12.7(1) 








Tavlor et al. f 1998) 


2451815.92 


2000.742 








12.294(23) 


12.256(31) 


12.189(23) 


2MASS 


2452568.53-2452568.67 


2002.805 








12.16-12.46 






Kato & Starkev (2002) 


2453651.87-2453651.95 


2005.770 






12.55(10)'^ 








Kafka et al. f2008) 


2453742.68-2454450.62 


2006.019-2007.959 








12.3(3)^^ 






AAVSO<^ 



"We have corrected the near-IR photometry listed in Table 2 of iHuber et al. I lll998l ) by using the 2MASS photometry of their comparison stars C(IR) and K(IR). 
These stars have 2MASS JHR^ magnitudes of 14.041(48), 13.772(62), 13.566(57) for C(IR) and 14.898(40), 14.304(58), 14.186(69) for K(IR), respe ctively. The 
2MAS S values are consi stently offset by 1. 0, 0. 3, and 1.1 mag in J, H, and Kb, respectively, compared to the comparison star photometry listed in iHuber et al.l 
jl998h (and also used bv lCiardi et al.lll998l and lTavlor et al.ll 19981 . 

^Mean V magnitude; the variations in the light curve range between V = 12.4-12.65. 

"^Photometry is "Visual", which is approximately /?-band. 

"^Henden, A. A., 2008, Observations from the AAVSO International Database, private communication. 



Table 3. Model component parameters 



Component Parameter Value Source 



System I 


Incliricition, i 


28+" 
^°-10 


Huber et al. (1998") 




Orbital Period, Porb (d) 


0.115063(1) 


Taylor et al. (1998) 




DistaincG, d (^pc^ 


364 


Taylor et al. (1998). this work 


WD: 


Tpmnpratnrp T j (K) 


45 000 


UlllO W Ui IV 




Mass M ^ (Mr-^) 


0.75 


this work 




Radius R j ( Rn,) 


0.0106 


this work 


SS; 


Sp6ctra,l Type 


M5.0 


Smith & Dhillon (1998'); Knieee (2006) 




Temperature, T'2 (K^ 


3030 


Cox (2000) 




Mass, M2 (Mq) 


0.21 


Cox (2000) 


ACD: 


Inner RaHins /? j ■ ( /? ^ 1 


1.0, 1.0 


this work 




Outer "Radius R ^ ^- ( R ^) 


35 35 


this work 




Critical Radius, i^crit (^wd) 


n/a, 12.25 


this work 




Inner Heiffht^ h ^ ( R ^) 


0.1, 0.1 


this work 




Outer Height'' , h^^ (R^d) 


0.1, 3.36 


this work 




Inner Temperature Profile Exponent, 7 


0.75, 0.75 


this work 




Outer Temperature Profile Exponent, 7out 


n/a, 0.70 


this work 




Mass Transfer Rate, M {Mq yr"!) 


1.1 X 10~«, 1.5 X 10-** 


this work 


CBD: 


Optical Depth Prescription 


thin 


this work 




Temperature Profile Exponent 


0.75 


this work 




Constant Height, /icbd (^wd) 


0.01 


this work 




Grain Density, Pgrain (g cm~^) 


3.0 


this work 




Grain Radius, rgjain (/^m) 


1 


this work 




Inner Radius, i?c:bd,in (^wd) 


700 


this work 




Outer Radius, Rcbd.out (^wd) 


15,000 


this work 




Inner Temperature, Tcbd in (K) 


500 


this work 




Outer Temperature, T^bd in (K) 


50 


this work 




Total Mass, M^bd (g) 


2.3 X 10^1 


this work 



00 
o 



''The disk height is the full thickness of the disk; that is, twice the distance from the disk midplane to its face. 



Note. — WD = white dwarf, SS = secondary star, ACD = accretion disk, CBD = circumbinary disk. Two values are listed 
for each of the ACD parameters, corresponding to the two models shown in Figure [2] all other parameters are the same for both 
models. 
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Fig. 1. — Per cent uncertainty (standard deviation) of measured IRAC flux densities as a 
function of average source flux density. The panels show (a) channel 1 (3.6 /im), (b) channel 
2 (4.5 /xm), (c) channel 3 (5.8 /xm), and (d) channel 4 (8.0 fim). The solid vertical line in 
each panel is the mean flux density of V592 Cas in the corresponding channel. The vertical 
dashed lines show the ±20% range over which uncertainty measurements were averaged to 
arrive at the scatter term uncertainties for V592 Cas (large diamonds). See §2.1.1l for details. 
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Wavelength (microns) 



Fig. 2. — Dereddened spectral energy distribution (top panel) and model (middle and bottom 
panels) for V592 Cas. The photometric data (see §3.ip are shown as filled circles. The system 
model (sohd line) is composed of a WD (dotted fine), M5.0 secondary star (short dashed 
line), steady state accretion disk (long dashed line), and circumbinary dust disk (dot-dot-dot 
dash line). The model in the middle panel utilizes a limb darkened fiat accretion disk with 
no irradiation, while the model in the bottom panel utilizes a limb-darkened fiared accretion 
disk with irradiation. See §3.21 for a discussion of the model parameters. 
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Fig. 3. — To-scale diagram of V592 Cas, including the model components (based on the 
parameters listed in Tableland discussed in the text). The diagram shows the secondary 
star (horizontal hatched area), WD (small filled circle), WD Roche lobe (dashed line), and 
the accretion disk (solid line around the WD marks the outer edge). The circumbinary dust 
disk, which has an inner radius of 700 -Rwd, is not shown. The inner and outer Lagrange 
points (plus symbols) and system center of mass (cross symbol) are also shown. The CV 
is depicted as viewed from "above" (i.e., with the orbital plane in the plane of the sky, 
equivalent to a system inclination of 0°). 



0.01 1 ... I ... I ... I ... I ... I ... I ... I ... I ... I ... I ... I ... I ... I . 
2 4 6 8 10 12 14 16 18 20 22 24 26 

Wavelength (microns) 

Fig. 4. — As in Figure [2], but showing a model that does not contain a circumbinary dust 
disk, in which a steady state accretion disk is the dominant component at all wavelengths. 
For clarity, only the total system model is plotted; the WD and secondary star are as shown 
in Figure O The top panel shows the effect of changing the outer radius of the accretion disk; 
from bottom to top the model curves have -Racd,out = 5, 15, 25, 35, 45-Rwd (the nominal value 
is 35-Rwd)- The bottom panel shows the effect of changing the mass transfer rate through 
the accretion disk; from bottom to top, the model curves have M = 0.01, 0.1, 1, 10, 100 times 
the nominal value of 1.5 x 10^^ M© yr~^. All other model parameters are as given in Table 

El 
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3 

Q 

o 




28 
26 
24 
22 
20 



1 1 1 1 1 1 1 


1 ' 1 ' 1 ' 1 ' 




(b) - 










I.I.I. 


1 . 1 . 1 . 1 , 



32 
28 
24 



3 4 5 

log(Time [yr]) 




, 1 , 1 1 1 


1 ' 1 ' 1 ' 1 ' . 




(d) : 

















Fig. 5. — Evolution of circumbinary dust mass as calculated using the toy model described 
in §4.61 In panels (a), (b), and (d), the CV's mass transfer rate evolves from lO"'''^ Mq yr~^ 
to 10^^ Mq yr~^ over 10^ yr. Panel (a) shows the effect of changing the dust dissipation 
efficiency, ei, with values of 0.01, 0.001, and 0.0001 (curves from bottom to top, respectively) 
while keeping the dust contribution efficiencies fixed at £2 = £3 = 0.001. Panel (b) shows 
the effect of changing the product of the dust formation efficiencies, 6263, with values of 
0.0001, 0.00001, and 0.000001 (curves from top to bottom, respectively) while keeping the 
dust dissipation efficiency fixed at ei = 0.001. In panels (c) and (d) the dust dissipation 
and formation efficiencies were fixed at ei = 62 = e3 = 0.001. Panel (c) shows the effect of 
changing the final mass transfer rate, with values of 10~^, 10~^, and 10^^^ Mq yr^^ (curves 
from top to bottom). Panel (d) shows the effect of changing the initial mass of circumbinary 
dust, with values of 10^^, lO^'', and 10^^ g (curves from bottom to top). Note the different 
vertical axis scale in panel (d). 



